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ABSTRACT 

>, 

O ! Recent photometric observations by the Hubble Space Telescope (HST) have 

^r^ \ revealed the physical properties of stellar galactic nuclei in nucleated dwarf galax- 

^ i ies in the Virgo cluster of galaxies. In order to elucidate the formation processes 

\Q \ of nucleated dwarfs, we numerically investigate gas dynamics, star formation, 

O 



and chemical evolution within the central 1 kpc of gas disks embedded within 
the galactic stellar components of non-nucleated dwarfs. We find that high den- 
_. , sity, compact stellar systems can be formed in the central regions of dwarfs as 

^ \ a result of dissipative, repeated merging of massive stellar and gaseous clumps 

j^ \ developed from nuclear gaseous spiral arms as a result of local gravitational in- 

stability. The central stellar components are found to have stellar masses which 
K^ \ are typically 5% of their host dwarfs and show very flattened shapes, rotational 

H \ kinematics, and central velocity dispersions significantly smaller than those of 

their host dwarfs. We also find that more massive dwarfs can develop more mas- 
sive, more metal-rich, and higher density stellar systems in their central regions, 
because star formation and chemical enrichment proceed more efficiently owing 
to the less dramatic suppression of star formation by supernovae feedback effects 
in more massive dwarfs. Based on these results, we suggest that gas-rich, non- 
nucleated dwarfs can be transformed into nucleated ones as a result of dissipative 
gas dynamics in their central regions. We discuss the origin of the observed cor- 
relations between physical properties of stellar galactic nuclei and those of their 
host galaxies. 

Subject headings: galaxies: dwarfs - galaxies: star clusters - globular clus- 
ters:general 
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Introduction 



A growing number of photometric studies of dwarf galaxies by the Hubble Space Tele- 
scope (HST) have recently revealed the physical properties of stellar galactic nuclei of dwarfs 
in the Virgo cluster of galaxies (e.g., De Propris et al. 2005; Grant et al. 2005; Ha§egan 
et al. 2005). i/^T photometric observations have also revealed that a significant fraction 
of late-type galaxies have very compact stellar systems in the very center of these galaxies 
(Phillips et al. 1996; Carollo et al. 1998; Matthews et al. 1999; Boker et al. 2004; Walcher 
et al. 2005). The physical properties of stellar galactic nuclei (e.g. masses and surface 
brightness) have been suggested to correlate with those of their host galaxies (Boker et al. 
2004). Spectrophotometric studies of nucleated dwarf galaxies have revealed the structural 
and kinematical properties of their central stellar components (e.g., Stiavelli et al. 2001; 
Geha et al. 2002; Lotz et al 2004). 

Although these recent observations, combined with previous ones on nucleated dwarfs 
(e.g., Binggeli & Cameron 1991; Ferguson & Binggeli 1994), have provided vital clues to 
understanding the formation of nucleated dwarf ellipticals (dE,Ns) and irregulars (dI,Ns), 
only a small number of theoretical works have so far investigated the formation processes of 
their stellar galactic nuclei (e.g., Tremaine et al. 1975; Oh & Lin 2000; Bekki et al. 2004). 
One of the most extensively discussed scenarios is that these stellar nuclei are formed from 
the merging of smaller stellar clusters of globular cluster-like mass in the central regions of 
galaxies (e.g.. Oh & Lin 2000; Bekki et al. 2004). These numerical works, which are based 
on dissipationless models of merging star clusters, have provided some predictions on the 
structural and kinematical properties of stellar nuclei. These previous works, however, did 
not discuss at all how gaseous dissipation and star formation are important in the formation 
of nucleated dwarfs, even though a number of spectrophotometric observations have already 
provided possible evidence for very young stellar populations in the stellar nuclei of some 
nucleated galaxies (e.g., Caldwell & Bothun 1987). 

The purpose of this Letter is thus to demonstrate, for the first time, that gaseous 
dissipation and star formation can play a vital role in the formation of nucleated dwarfs (i.e., 
dE,Ns and dI,Ns). Based on the results of chemo dynamical simulations of the gas dynamics 
and star formation within the central Ikpc of non- nucleated dwarfs, we demonstrate that non- 
nucleated dwarfs can be transformed into nucleated ones as a result of dissipative nucleation. 
We investigate the structural, kinematical, and chemical properties of stellar galactic nuclei in 
nucleated dwarfs with different masses and gas mass fractions and thereby provide theoretical 
predictions that can be tested against recent observations. Although the present study 
suggests the importance of gas clump formation in inward mass-transfer processes for dwarfs, 
this importance has already been proposed by Shlosman & Noguchi (1993) and Noguchi 
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(1999) in the context of triggering mechanisms of starbursts and bulge formation. 



2. The model 

We investigate the dynamical evolution of stellar and gaseous components within the 
central 1 kpc of a galaxy, via chemodynamical simulations carried out on a GRAPE board 
(Sugimoto et al. 1990). Since the numerical methods and techniques we use for the GRAPE 
system have already been described in detail elsewhere (Bekki & Shioya 1999), we give only 
a brief review here. 

Old stellar components of dwarf galaxies (dis and dEs) are observed to have exponential 
luminosity profiles (e.g., Ichikawa, Wakamatsu, & Okamura 1986) and therefore the projected 
radial density profile of the stellar component with mass Mgph and scale length Ogph is as- 
sumed to be an exponential one with the central surface density of Sgph^o- Guided by recent 
observations on the scaling relation between Ssph,o and Mgph (i.e., Ssph,o oc Mgph^'^^; Kauff- 
mann et al 2003), we adopt the Mgph — agph relation of agph = CqX Msph°'^^, where Cq is chosen 
such that the mass and the scale length of the Galactic exponential bulge can be reproduced 
reasonably well. For example, this equation gives Csph = 200pc for Mgph = lO^M©. 

The gas disk in a dwarf is modeled as a collection of discrete gas clouds with the observed 
mass-size relationship of GMCs (Larson 1981). The thin gas disk composed of these clouds 
is assumed to have an exponential profile with the scale length being the same as Ogph- The 
mass ratio of the gas disk to the initial stellar component is regarded as a free parameter (/g) 
that ranges from 0.02 to 0.5. Every pair of two overlapping gas clouds is made to collide with 
the same restitution coefficient (Hausman & Roberts 1984) and lose its random kinematic 
energy through gaseous dissipation. 

Field star formation is modeled by converting the collisional gas particles into collision- 
less new stellar particles according to the algorithm of star formation described below. We 
adopt the Schmidt law (Schmidt 1959) with exponent 7 = 1.5 (1.0 < 7 < 2.0, Kennicutt 
1998) as the controlling parameter of the rate of star formation. These stars formed from gas 
are called "new stars" (or "young stars" ) whereas stars initially within a spheroid are called 
"old stars" throughout this paper. Chemical enrichment through star formation and super- 
novae feedback during the stellar nucleus formation and evolution is assumed to proceed 
locally and the values of the return parameter and the chemical yield are set to be 0.3 and 
0.005, respectively. About 10% of the total energy of one type-II supernovae (~ 10^^ erg) is 
assumed to be converted into kinematic energy of gas around the supernovae (Thornton et 
al. 1998) for the Salpeter IMF. The feedback effects are stronger in low-mass dwarfs for the 
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adopted mass-size relation in the present study. 

We mainly investigate projected radial density profiles (/Xg) of new and old stars in the 
central regions of spheroids. In order to compare the simulated profiles (/ig) with the observed 
i?— band surface brightness profiles (/iB), we use the formula of /ib = 27.05 — 2.51og(/B), 
where Jb is the 5— band surface brightness of a stellar system measured in units of L0 pc^^. 
We estimate Jb from /ig for a given mass-to-light-ratio {M/ Lb) in the B— band. We assume 
that M/Lb for old stars [{M/Lb)^-^^ is 10, corresponding to a single stellar population 
(SSP) with age 10 Gyr and solar metallicity (Vazdekis et al. 1996), and that for new stars, 
(M/Lb)^^^ is 0.5, corresponding to an SSP with age 0.5 Gyr and solar metallicity. 

We firstly describe the results of the fiducial model with Mgph = IO^Mq and /g = 
0.2, which shows typical behaviors of dissipative nucleus formation in dwarfs. Then we 
describe briefly the dependences of the results on Mgph and /g. The details of the parameter 
dependences (including those on M/ Lb) will be given in our forthcoming papers. The total 
number of particles is 10^ for the old stars (A^ow) and 2 x 10"^ for the initial gas (iVgas), and 
the time integration of the equation of motion is performed by using the 2nd-order leap-frog 
method. 



3. Result 

Figure 1 shows how a compact stellar system is formed in the central region of the 
star-forming gas disk embedded in the old stellar spheroid for the fiducial model. Numerous 
spiral arms composed only of gas clouds first develop owing to the stronger self-gravity of 
gas (/g = 0.2) as the gas disk rotates (T = 62 Myr). Because of the enhanced cloud-cloud 
collision rate and the resultant efficient gaseous dissipation within the gaseous spiral arms, 
small gas clumps composed of several tens of gas clouds can be formed within each spiral arm 
(T = 250 Myr). Some of these small gaseous clumps are massive, compact, and strongly self- 
gravitating so that gas can be continuously converted into stars within the clumps (T = 250 
Myr). These clumps can spiral in owing to dynamical friction against the stellar background 
of old stars (i.e., the stellar spheroid), and consequently merge with one another to form two 
very massive clumps composed mostly of new stars (T = 250 Myr). 

The two massive clumps that develop within the central 200 pc can interact with each 
other to form tidal arms composed only of new stars (T = 374 Myr). Their mutual tidal 
interaction results in the loss of orbital angular momentum, and they therefore finally merge 
with each other to form a single compact stellar system in the central region of the stellar 
spheroid. Although the single nucleus is initially not in the very center of the spheroid. 



-5- 



dynamical friction between the nucleus and the stellar background can finally transfer the 
nucleus into the mass center of the spheroid within ~ 100 Myr. Thus a stellar galactic 
nucleus composed of new stars can be formed from merging between massive stellar and 
gaseous clumps in the central region of the spheroid. 

Figure 2 clearly shows that the non-nucleated dwarf is transformed into a nucleated 
system with a very high central surface brightness owing to dissipative formation of the 
compact stellar system in its center. The derived large differences in /x^ between old and new 
stars results from the combined effect of the central high stellar density (due to dissipative 
stellar nucleus formation) and the small M/ L of the young stellar population in the nucleus. 
Figure 3 shows the radial profiles of the velocity dispersion (a) estimated separately for the 
old stars and the new stars at T = 686Myr. It is remarkable that the new stars show a 
velocity dispersion profile that increases radially from a central value of a^ ~ 15 km s^^, 
whereas the old stars show a nearly constant profile with a^ of ~ 60 km s^^. The essential 
reason for the lower o"o in the new stars is that the stellar nucleus, which is formed from 
dissipative merging of clumps, can be strongly self-gravitating so that ctq is determined by 
the nucleus itself (i.e., by its smaller mass of ~ 0.05 in our simulation units). The stellar 
nucleus shows a significant amount of rotation (20 — 50 km s~^ within the central lOOpc) 
and a very fiattened intrinsic shape. 

Figure 4 shows the age-metallicity relation of the new stars that are located within the 
central 50 pc (thus within the nucleus) of the galaxy at T = 686Myr. The mean age and 
metallicity of this stellar nucleus are estimated to be 3.7 x 10^(= 10^'^^) yr and —0.29 in 
[Fe/H], respectively. The final metallicity of —0.29 is significantly higher than the initial value 
([Fe/H]i = —0.44) for the chemical yield of 0.005 and [Fe/H]i = —0.44, which means that 
the stellar nucleus is more metal-rich and younger than the background stellar (spheroidal) 
components. Figure 4 clearly indicates that the stellar nucleus is composed of stars with 
different ages and metallicity: The ages of the new stars ranges from to 6.3 x 10^(= 10^'^) 
yr to 6.8 x 10^(= 10^'^) yr whereas the metallicities range from —0.44 to —0.077 in [Fe/H]. 
This wider range of ages and metallicities reflect the fact that the stellar nucleus can form 
from merging between different clumps, which have different star formation and chemical 
evolution histories. It is clear from Figure 4 that the more metal-rich stars are more likely 
to be younger. 

The models with different M^ph (2.5 x lO^M©, 1.0 x lO^M©, 5.0 x lO^M©, and 1.0 x 
lO^M©) and different /g (=0.02, 0.05, 0.2, and 0.5) are investigated in the present study. 
The dependences of the physical properties of the compact stellar systems (which can be 
regarded as stellar nuclei) on Mgph and /g are briefly summarized as follows: Firstly, the 
total masses of the stellar nuclei (Mnuc) in units of M© within 50 pc (5% of an initial gas 
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disk size) is correlated with the initial stellar masses (Mgph), and this correlation can be 
described as M^uc = 0.046 x Mgph — 0.0030. Here we apply a least-squares fit to the results 
of the models that have different Msph and clearly show stellar nucleus formation. If we 
make a least-squares fit to the logj^gMnuc — logio^sph relation, we can derive the relation 
of logigMnuc = 1-1 X logioMsph — 2.5. Secondly, the mass fraction of the stellar nucleus 
(/nuc) in nucleated dwarfs is correlated with Mgph and the correlation can be described as 
/nuc = 0.021 X (yq9^) +0.023. Thirdly, the metallicities of stellar nuclei are higher than those 
of their host galaxies in most models. Accordingly stellar nuclei can be significantly redder 
than their host galaxies after significant aging of their stellar populations. More massive 
dwarfs can have more metal-rich stellar nuclei. 

Fourthly, stellar nuclei do not develop within a few Gyr in the models with smaller /g 
(< 0.05), which strongly suggests that self-gravity of gas is quite an important factor for 
dissipative formation of stellar nuclei. The central surface brightness (/iq) of stellar nuclei 
is higher in models with higher /g, which implies that the observed diversity in ^q between 
different dwarfs (e.g., De Propris et al. 2005) can be understood in terms of a different 
amount of gas within the dwarfs. Fifthly, dense stellar nuclei are more likely to develop 
in dwarfs with higher stellar densities owing to their deeper gravitational potentials (i.e., 
weaker supernovae feedback effects). 



4. Discussions and conclusions 

It has long been known that more luminous dwarf galaxies (dE) are more likely to contain 
stellar nuclei (e.g., van den Bergh 1986): The number fraction of dE,Ns among dEs is almost 
100 % for Mb ~ — 18mag and about 20 % for Mb ~ — 12mag. The present simulations 
have demonstrated that less luminous galaxies can have less remarkable stellar nuclei and 
some stellar nuclei can hardly be identified observationally as distinct nuclei. They also have 
suggested that the origin of this dependence is closely associated with the fact that stellar 
nucleus formation is more strongly suppressed by stronger feedback effects in less luminous 
galaxies. These simulation results suggest that the origin of the observed dependence of 
dE,N fraction on the host luminosities can be understood in terms of the luminosity (or 
mass) dependence of the effectiveness of supernovae feedback in stellar nucleus formation. 

Recently Lotz et al. (2004) have revealed that some stellar nuclei of dE,Ns have colors 
that are significantly bluer than those of their host galaxies and can be due to recent star 
formation episodes. Several previous spectroscopic observations of stellar nuclei in dE,Ns 
also revealed possible evidences for younger stellar populations these stellar nuclei (e.g., 
Caldwell & Bothun 1987; Bothun & Mould 1988). The present study suggests that if these 
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apparently young stellar nuclei were formed dissipatively from gas, as our simulations have 
demonstrated, they should be significantly more metal-rich than their host galaxies. It is 
thus doubtlessly worthwhile for future spectroscopic observation to investigate the metallicity 
differences between stellar nuclei and their host galaxies. 

We have demonstrated that gas dynamics and star formation in the central 1 kpc of 
dwarfs can be important for the transformation of non-nucleated dwarfs into nucleated sys- 
tems. One of advantages of the dissipative nucleus formation scenario is that it can naturally 
explain how massive stellar and gaseous clumps, which are "building blocks" of stellar nu- 
clei, can be formed in the central regions of gas-rich dwarfs. This scenario, however, has not 
yet provided clear explanations on some of the key observed correlations between dynamical 
properties of stellar nuclei and those of their host galaxies (e.g., Boker et al. 2004). 

We are grateful to the anonymous referee for valuable comments, which contribute to 
improve the present paper. K.B. and W. J. C acknowledge the financial support of the 
Australian Research Council throughout the course of this work. The numerical simulations 
reported here were carried out on GRAPE systems kindly made available by the Astronomi- 
cal Data Analysis Center (AD AC) at National Astronomical Observatory of Japan (NAOJ). 
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Fig. 1. — Morphological evolution of gas (cyan) and new stars (magenta) in the fiducial 
model projected onto the x-y plane. The time T shown in the upper left corner of each panel 
represents the time that has elapsed since the simulation starts. One frame measures 2 kpc 
for the upper four panels and 400 pc for the lower ones. 
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Fig. 2. — 5-band surface brightness profiles (/xb) for old stars (dotted), new stars (thin 
solid) and total (i.e., old + new stars, thick solid) at T = 686Myr in the fiducial model 
with (M/Lb)oi<j = 10 and (M/Lb)^^^ = 0.5. The results correspond to the epoch when the 
simulated nucleated dwarf is young (~ 0.5 Gyr). 
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Fig. 3. — Radial profiles of velocity dispersion for old stars (thin solid) and new ones (thick 
solid) at T = 686Myr for the edge-on view (i.e., projected onto the x-z plane) in the fiducial 
model. 
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Fig. 4. — Distribution of new stars on the age-metallicity map at T = 686Myr in the fiducial 
model with the initial gaseous metallicity [Fe/H]j of —0.44. Here only new stars within the 
central 50 pc (thus those regarded as being in the stellar nucleus) are plotted. The horizontal 
and vertical dotted lines represent the mean metallicity and age, respectively. 



